Astronomy 8l Astrophysics manuscript no. coalsack'071 1 
July 13,2011 



©ESO2011 



The Coalsack near and far 

H. Beuther\ J. Kainulainen\ Th. Henning\ R. Plume^, and F. Heitsch^ 

^ Max-Planck-Institute for Astronomy, Konigstuhl 17, 69117 Heidelberg, Germany, e-mail: beuther@mpia.de 

^ Department of Physics and Astronomy, University of Calgary, Calgary, Canada 

^ Department of Physics & Astronomy, University of North Carolina at Chapel Hill, CB 3255, Chapel Hill, NC 27599-3255, USA 

Preprint online version: July 13, 2011 



o 

1—5 

<N 



Oh' 

O 



> 

00 
On 
(N 
(N 

• 

O 



X 



ABSTRACT 

Context. The large Coalsack dark cloud is one of the most prominent southern starless clouds, which is even visible to the naked eye. 

Furthermore, it is one of the rare molecular clouds without clear signs of star formation. 

Aims. We investigate the dynamical properties of the gas within the Coalsack. 

Methods. The two highest extinction regions were mapped with the APEX telescope in ^^C0(2-l) comprising a region of ~1 square 

degree. 

Results. In addition to the well-known, nearby gas component around -4 km s~\ we identified additional molecular gas components - 

in particular a second extended molecular cloud at a velocity of ~ -30 km s~^ and an estimated distance of ~3. 1 kpc - that dominate the 

column density and visual extinction distributions in the northeastern part of the Coalsack. Although comprising ~2600Mo, the mass 

of this distant cloud is distributed over an extent of ~73 pc, much larger than typical high-mass infrared dark clouds. Its filamentary 

structure is consistent with a compressible gaseous self-gravitating cylinder, and its low mass per length indicates that it may be 

stable against gravitational collapse. We find barely any mid-infrared emission in archival MSX data, which is indicative of almost no 

star-formation activity in the near and far cloud complexes. The nearby clouds have narrow, almost thermal velocity dispersions with 

median values between 0.2 and 0.4kms~\ which is also consistent with low star-formation activity. Only Tapia's Globule 2 exhibits 

a velocity dispersion increase toward the extinction peak and peak- velocity gradients over the core, which is indicative of a state of 

elevated dynamical properties. 

Conclusions. The Coalsack is not one single coherent structure, but consists of several cloud complexes nearby as well as at several 

kpc distance. All studied clouds appear as starless low-turbulence regions that may not even collapse in the future. Only one globule 

exhibits more dynamical signatures and is a good candidate for present/future star formation. 

Key words. Stars: formation - ISM: individual: Coalsack - ISM: clouds - ISM: kinematics and dynamics - ISM: structure 



1. Introduction 

■ The Coalsack is one of the closest and most widely-recognized 
dark clouds in the southern hemisphere (between Galactic lon- 
gitudes 3 05 > / > 300 deg) because it i s visib le to the bare 
eye (e.g.jNymanl l2008l) . lKainulainen et alJ (l2009h recently pro- 

, duced extinction maps based on 2MASS near-infrared data of 
all supposedly nearby molecular clouds including the Coalsack. 
Figure[T] shows the large-scale extinction (A^) map of this region. 
While the visual extinction through the Coalsack is generally 
low, there are a few regions where the extinction A^ rises above 

I 6 and even 10 mag. Although the extinction maps set significant 
constraints on the morphology and column density distributions, 
allowing us to derive results about cloud structure and stabil- 
ity a s well as potential column de nsity star-formation thresh- 
olds ("Kain ulainen et al.l 120091 12011) , these data lack any kine- 
matical information about the cloud. Furthermore, this cloud is 
one of the rare cases of a nearby molecular cloud de void of any 
star-f ormation signatures (e.g., Nyman.2008,: Kainulainen et al.l 
l2009l) . This is one of the reasons why it has not received as much 
attention in the past as typical star- forming clouds like Taurus of 
p Ophiuchus. Despite this, the Coalsack therefore constitutes an 
excellent laboratory to study pristine and undisturbed molecular 
clouds. 
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Earlier s pectral-line rnapping on the Coalsack was conducted 
in ^^CO bv (Nvman et al.1 (Il989l) at a spatial resolution of 8.8', 
and in ^^CQ bv lKato et al.l (Il999h at a resolution of 2.1'. Thus, 



Fig. 1. 2MASS extincti on map of the whole Coalsack by 
iKainulainen et al.l (l2009l) . The contour levels are at 2, 6 and 
10 mag extinction. The two boxes outline the two regions 
mapped in i^C0(2-l) and C^^0(2-l) for this work. 



even basic dynamical measures (e.g., central velocities, dis- 
persions, etc.) were limited to spatial scales larger than what 
is expected for molecular clumps or cores (or equally limited 
in dynamical range/velocity resolution). The typical velocities 
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of the region are between -6 and Okms ^ (e.g., iNyman et al.l 
Il989l) . and the distance estimates vary between 100 and 200 pc 
(e.g., [Franco 1995; Corrad i etal. 1997; Knude & Hog 1998; 
lNvmanlf2008h . iNvmanetaD (Il989i) reported the presence of an 
additional velocity component at -35kms"\ which they at- 
tribute to a background cloud associated with the Sagittarius- 
Carina arm that is unlikely to be responsible for any of the 
extinction features. Some individual fields of the Coalsack 
have been mapped in tracers of higher density by some au- 
thors, mainly concentrating on the region known as "Tapia's 
Globule 2", whi ch presumably is one of the few dense cores in 
the region (e.g., Tapia 1973; Ladaetal. 2004; Rathborne et al. 
[2009). Furthermore. Font ani et al. (2005), Beltran et al. (2006) 
and iMiettinen & Harjul (12010) studied a smaller sub-filament 
just adjacent to the eastern edge of our northeastern field. These 
studies, however, cover only a minuscule portion of the complex 
and therefore are not very informative for the dynamics of the 
complex itself. 



T.b[K] 
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Fig. 2. Average ^^C0(2-l) spectra toward the northeastern re- 
gion Rl (gray with this histogram) and the southwestern region 
R2 (thick histogram multiplied by 0.2 to be on the same scale). 

To fill this missing link was one of the motivations mapping 
parts of the Coalsack in spectral line emission. For this pur- 
pose, we selected the two major regions with significant frac- 
tions of cloud structure above ~6 mag extinction, marked by the 
two boxes in Figure [T] The southwestern region also comprises 
"Tapia's Globule 2". These two regions cover in total about 1 
square degree and were mapped with APEX in ^^C0(2-l) and 
C^^0(2-l) at a spatial resolution of ~ 27.5''. 

2. Observations and data 

The Ci^0(2-1) and ^^C0(2-l) data at 219.560GHz and 
220.399 GHz were observed sim ultaneously with the Atacama 
Pathfinder Experiment (APE:x[3, iGtisten et all l2006h between 
April and June 2010 in the 1 mm band for the two regions Rl 
and R2 marked in Fig. [T] in on-the-fly mode. The APEXl re- 
ceiver of the SHeFI receiver family has rec eiver temperatures 
of -130 K at 220 GHz ( Vassilev et al.1. [2008) . and the average 
system temperatures during the observations were -220 K. Two 
Fast-Fourier-Transform- Spectrometer (EFTS, iKlein et ani2006h 
were connected, covering ~2GHz bandwidth between 219 and 
221 GHz with a spectral resolution of ~0.17kms"^ The data 



^ This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 



were converted to main-beam brightness temperatures Tmb with 
forward and beam efficiencies at 220 GHz of 0.97 and 0.82, 
respectively (Vassilev et al., 2008). The average Icr rms in the 
maps of emission-free channels is ~0.75 K. Moment maps were 
conducted clipping all pixels below a 2 K threshold. The selected 
OFF-positions were at the edge of the Coalsack and were free of 
apparent CO emission (R.A. (J2000) 12:57:20.514 Dec. (J2000) 
-60:59:43.42 for Rl and R.A. (J2000) 12:31:51.0 Dec. (J2000) 
-63:57:56.9 for R2). The FWHM of APEX at the given frequen- 
cies is ~ 27.5''. 

The near-inf ra red extinction map is taken from 
Kainulai nen et al.l (l2009l) . The spatial resolution of this 
map is ~ 90''. 



3. Results 

3.1. The Coalsack: Near and far components 

The ^^C0(2-l) emission in the northeastern region Rl ex- 
hibits several velocity components (see average spectrum in 
Fig. [2]), most prominently the well-known Coalsack component 

at 2kms"^ as well as a second even stronger component at 

30kms"^ Furthermore, there are additional components at 

more negative velocities, one centered at 38kms"^ and one 

at 44kms"^ The south-eastern region R2 exhibits only one 

velocity component at 4 km s"^ . The little velocity off'sets be- 
tween the 2kms"^ and 4kms"^ components of regions 

Rl and R2 imply a small large-scale velocity gradient over the 
entire Coalsack complex. How do these components now corre- 
spond to the extinction in the Coalsack? 

Figures [3] and [4] present integrated ^^C0(2-l) maps of the 
strongest emission features toward the northeastern region Rl 
and the southwestern region R2, respectively. While the ~ 
-4kms"^ component in region R2 exhibits the expected struc- 
ture, peaking in th e direction of the extinction peak of Tapia's 
globule 2 (see also lLada et"aDl2004l: [Rathborne et al.|[2b09.) . the 
region Rl clearly exhibits two spatially very distinct veloc- 
ity components around ~ -2 and ~ -30kms"^. Furthermore, 
Figure [5] presents the integrated ^^C0(2-l) maps of the addi- 
tional weaker components at ~ -38kms"^ and ~ -44kms"^. 
While we can identify spatial structures from these two velocity 
components as well, they are less prominent and also less ex- 
tended. Therefore, we omit the latter two velocity components 
in the following discussion. 

The C^^ 0(2-1) emission is comparably weak and we only 
barely detect some structures. While we detect C^^ 0(2-1) in 
the ~ -4kms"^ component toward R2, we only detect the 
~ -30kms"^ component toward Rl. Therefore, here we only 
work with th e^^C0(2-l) d| ata. 

Although iNyman et al.l (1 19891) reported spectra with a sec- 
ond velocity component, their mentioned positions were far off'- 
set from the northeastern region Rl. In addition to this. Figure 
[3] clearly shows that the extinction is not caused by the nearby 
component at the velocity of ~ -2kms"^ but is mainly as- 
sociated with the second much stronger velocity component at 
30kms"^ This is important because it implies that the vi- 
sual dark cloud Coalsack is not one coherent low-mass cloud as 
believed so far, but that it is comprised of several clouds at dif- 
ferent distances, all combined to produce the dark features on 
the sky. The dark features come primarily from the distant com- 
ponent. 

Following the rot ation curves by either lBrand & Blita (Il993h 
or iReid et al.l (l2009h . for the -30kms"^ component the kine- 
matic near distances are 2.7 or 3.1 kpc, respectively. Because 
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Fig. 3. Compilation of the integrated emission in grayscale (with thin contours) from the -2kms"^ and -30kms"^ components in 
^^C Q(2-1) toward the northea stern region Rl in the top and bottom panels, respectively. The thick contours outline the extinction 
from lKainulainen et alJ (l2009l) in 1 mag steps between Ay of 5 and 9 mags plus additional contours of 1 1 and 1 3 mags. The integration 
regimes for the top and bottom panels are [-4,2] and [-35,-24] km s"^ respectively. The thin contours for the -4kms"^ map start at 

1 Kkm s"^ and continue in 1 Kkm s"^ steps, whereas for the -30 km s"^ component the contours start at 2 Kkm s"^ and continue in 

2 Kkm s"^ steps. A linear scale-bar is shown in both panels 



we see the Coalsack in extinction, we can safely take the near 
distance since at far distances one would not see any extinc- 
tion. For comparison jFontani et al. ' (^2005) derived a distance of 
2.44 kpc for the position of the IRAS 13039-6108 source from 
CS data. They used a velocity of -26.2kms"^ consistent with 
our data at the eastern edge of the map (see Fig. [7]), as well as the 
iBrand & Blitzl (|T993) rotation curve. Because we use an average 
velocity over the whole structure for the distance determination 
(Fig. [7]) and try out different rotation curves, we will stick to the 
distance of 3.1 kp c derived with the most recent rotation curve 
(iReid et al.L 12009 ). Qualitatively speaking, these results imply 
that not only a small sub-filament around the IRAS source, but 
the whole northeastern high extinction region Rl is not, as usu- 
ally assumed, a low-mass dark cloud at distances below 200 pc, 
but is a more massive molecular cloud at a distance of several 
kpc. 



Employing the Galactic structure outlined in iReid et afl 
(12009), our lines-of- sight pass through the near and far sides 
of the Galactic spiral arms Carina- Sagittarius and Crux-Scutum 
(<2pc and >3kpc distance, respectively). While the Carina- 



Sagittarius arm should accordingly be responsible for the back- 
ground structures of the nearby cloud components, likely the 
Crux-Scutum arm - in which the distant cloud component 
should be embedded - also provides the emission against which 
to produce the observed extinction features. 

Assuming optically thin ^^C0(2-l) line emission, we can 
calculate the gas masses for the diff'erent components follow- 
ing Rohlfs & Wilson (2006). The -30km s"^ component in re- 
gion Rl is assumed to be at a distance of 3.1 kpc (labeled Rlfar 
from now on), whereas for the -2/-4kms"^ components in Rl 
and R2 we use a distance of 175 pc (and label it Rlnear in 
the northeastern region Rl). Because all clouds are extinction 
features and hence are largely comprised of cold, mainly star- 
less gas, an average temperature of 15K appears reasonable. 
With these numbers, we derive a mass of ~2600Mo for the - 
30kms"^ Rl cloud, whereas the two nearby -2/-4kms"^ of Rl 
and R2 only have masses of -1.2 and ~1O.4M0, respectively. 
Hence the 30kms"^ cloud at a distance of 3.1 kpc is 2-3 orders 
of magnitude more massive than the other two components at 
the more typical Coalsack distance. However, one should keep 
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Fig. 7. Compilation of ^^C0(2-l) 1st moment maps toward the northeastern region Rl as marke d in the panels. Th e top an d bottom 
panel show the -2 and -30kms"^ respectively. The thick contours outline the extinction from iKainulainen et al.l (l2009l) in 1 mag 
steps starting between Ay of 5 and 9 mags plus additional contours at 1 1 and 13 mags. Linear scale-bars are shown in the bottom-left 
corner of each panel. 
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Fig. 4. The grayscale with thin contours shows the in- 
tegrated ^^C0(2-l) emission toward the southwestern re- 
gion R2. Contour levels start at 2Kkms"^ and continue in 

1 Kkms"^ steps. The thic k contours outline the extinction from 
IKainulainen et al.l (l2009l) in 1 mag steps starting at an Ay of 

2 mags. A linear scale-bar is shown in the bottom-right corner. 



in mind that a -2600 M© cloud, that is distributed over a size of 
~ 73 pc in east- west direction, cannot be considered as a typical 
high-mass star-forming region. The mass is quite smoothly dis- 
tributed over a large spatial area. For comparison, cloud masses 
of typical low- and high-mass star-forming regions like Taurus 
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Fig. 5. Integrated ^^C0(2-l) images in the northern Rl region of 
the two other velocity components seen in Fig.[2j The grayscale 
shows the integration regime [-41,-36] km s"^ and the contours 
present the velocity regime [-50,-41] km s"^ 



or Orion vary between s everal thousand to several 10^ M© (e.g., 
IKainulainen et al.ll2009l) . 

To check for signatures of ongoing star formation, we in- 
vestig ated the 8 and 21yum images of the MSX satellite dPricel 
1 19951) . PI. Figure [6] presents overlays of the 8yum emission with 
the extinction map. The northern region Rl shows one bright 
8//m source with a bit of extended diffuse emission in the 
west. This bright region may be associated with the extinction 



^ Our Coalsack targets are only barely covered by Spitzer because 
they are at relatively high and low galactic longitude, and GLIMPSE 
only covered latitudes < +1 deg dChurchwell et al.i.i2009i) 
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Fig. 6. The grayscaling presents the 8//m data toward the Rl 
(top) and R2 (bottom) regions as observed with MSX The con- 
tours outline the extinction from iKainulainen et alJ (I2QQ9 ) for 
Rl in 1 mag steps between Ay of 5 and 9 mags plus additional 
contours at 1 1 and 13 mags, and for R2 in 1 mag steps starting at 
an Ay of 2 mags. The white stars mark the only sources that were 
also detected with MSX at 21 yum. 



gradient implies that this is one coherent structure, even the parts 
with lower extinction and almost no ^^C0(2-l) detection in the 
middle of the filamentary structure. 
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Fig. 8. Outline of the clump distribution for R lfar. The grayscale 
with contours shows the extinction map from IKainulainen et al.l 
(2009) in 1 mag steps between Ay of 6 to 9 mags plus addi- 
tional contours at 11 and 13 mags. A linear scale-bar is shown 
at the bottom-left. The east- west extent of this extinction map 
is slightly larger than that of the spectral line data. The trian- 
gles mark the c lump peak positions as derived via gaussclump 
dWilliams et al.L [1994), and the connecting line visualizes the 
separations measured between the clumps. The histogram above 
shows the distribution of nearest neighbor clump separations as 
derived from the extinction map. 



peak at an approximate separation of 2.6'', however, no 21yum 
emission is detected from that region. While the 8yum peak is 
neither clearly associated with the -30kms"^ nor with the - 
2kms"^ ^^CO component, the more diff'use 8yum emission ex- 
tending from this peak toward the north-west is associated with 
near-infrared extinction as well as ^^C0(2-l) emission of the - 
30kms"^ component. Therefore, it appears more likely that the 
Syum peak is also associated with the -30kms"^ component. 
In general, there is barely any 21yum emission detected toward 
both regions, implying no or only a marginal level of embedded 
star formation. This is consistent with the general picture of the 
Coalsack being largely a starless cloud. 



3.2. Kinematic analysis of the Coalsack clouds 

In addition to the distance and mass estimates, the ^^C0(2-l) 
data also allow us to characterize the kinematic properties of the 
two - or more correctly three - regions. Here, we will focus on 
the northeastern regions (Rlfar and Rlnear) and the southwest- 
ern region R2 separately. 

3.2.1 . The more massive and more distant cloud Rlfar 

Figure [7] presents the ^^C0(2-l) first-moment maps (intensity- 
weighted peak velocities) of the two northeastern regions Rlfar 
and Rlnear. Covering an east- west length of ~ ST (or ~ 
1.35 deg) equivalent to a linear extent -73 pc, the velocity gra- 
dient across this structure is relatively small with a velocity dif- 
ference of only ~8kms"^ (from 26 to 34kms"^), cor- 
responding to ~0.1kms"^ pc"^ The relatively smooth velocity 



Inspired by the recent work of I Jackson et al.l (l201Qh on 
the large (~80pc length) "Nessie" filamentary infrared dark 
cloud (IRDC), we investigated whether the -30kms"^ com- 
ponent Rlfar may also be consistent with a gravitationally 
bound gaseous cylinder that may fragment because of fluid in- 
stabilities. For that pu rpose we used the clumpfind algorithm 
(IWilliams et al.L Il994l) to extract clump peak positions from the 
near-infrared extinction map. Because we are only interested in 
the peak positions here, we started with relatively high 6 mag 
contours and continued at 1 mag levels as input parameters for 
clumpfind. Because the filamentary structure is a bit larger than 
our ^^CO coverage, we used also a slightly larger extinction map 
for the purpose of the filament analysis (Fig. (8]). The projected 
east- west extent of the analyzed structure is -88 pc. The ex- 
tracted positions are shown in Figure [H Based on these peak 
positions we calculated the projected linear separations for these 
peaks assuming a distance of 3.1 kpc. One should keep in mind 
that all observed separations discussed below have always to be 
considered as lower limits because we have no information about 
the inclination angle of the filament with the plane of the sky. 
Figure [8] also presents a histogram of the distribution of linear 
separations, and the peak of the distribution is at about ~3.5 pc. 
While the mean separation is ~4.6pc, the median is at ~4.1 pc. 
This separation is resolved well by the spatial resolution of the 
extinction map of ~ 1 .35 pc, and random placement of the cores 
within the filament would produce a very diff'erent distribution 
with a mean core separation of ~2pc. The median separation 
within the Rlfar fil ament is consi s tent w ith that derived for the 
"Nessie" IRDC by iJackson et al.l (l2010l) . For more discussion 
see section |4J1 
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Fig. 9. Compilation of ^^C0(2-l) 2nd moment maps toward the northeastern region Rl as mark ed in the panels. Th e top and bottom 
panel show the -2 and -30kms"^ respectively. The thick contours outline the extinction from iKainulainen et al.l (|2009) in 1 mag 
steps between Ay of 5 and 9 mags plus additional contours at 1 1 and 13 mags. Linear scale-bars are shown in the bottom-left corner 
of each panel. 



Figure [9] presents the 2nd moment maps - or the intensity- 
weighted velocity dispersion. While there is a spread in velocity 
dispersion between approximately 0.4 and 4 km s"^ over the en- 
tire region, most observed line-widths are at the lower end of 
the distribution with an average velocity dispersion over the en- 
tire field of ~1.26kms"^ and a median velocity dispersion of 
~0.8 km s"^ . This relatively narrow average line- width compared 
to more evolved star-forming regions confirms the previous as- 
sessment from the near- and mid-infrared data, that the whole 
region exhibits low star- formation activity capable of contribut- 
ing to the line- width broadening. 



3.2.2. The low-mass nearby cloud Rinear 

Similar to Rlfar, Figures [7] and [9] also present the 1st and 2nd 
moment maps of the -2kms"^ component in the Rl region. 
From now on we call that Rinear. As already outlined in sec- 
tion [Til there is barely any spatial association of the extinction 
map with the -2 km s"^ ^^C0(2-l) emission. The velocity spread 
of the region between R.A. 13:00:00 and the western edge of 
the map is only about 4kms"^ between ~ -3 and ~ Ikms. 
This velocity gradient covers an east- west extent of ~ 32', cor- 
responding at a distance of 175 pc to an approximate linear ex- 
tent of ~ 1.6pc. This results in an approximate velocity gradi- 
ent of 2.6 km s"^ pc"^ While this velocity spread appears to be 
relatively small compared with the linear velocity gradient one 
can derive for Rlfar (~0.1 kms"^ pc~^), it is still comparatively 



large. Or the other way round, the velocity gradient for Rlfar is 
extremely small over a much larger spatial extent (see previous 
section). 

The 2nd moment or velocity dispersion distribution for 
Rinear is even narrower than that for Rlfar (Fig.©. While there 
are a few small regions exhibiting velocities >1 km s"^ the aver- 
age and median velocity dispersions for Rinear are 0.85 km s"^ 
and 0.2kms"\ respectively. While the average represents the 
spread in velocities a bit better, the median clearly shows that 
most observed line-widths are very narrow. At an approximate 
cold temperature of 15K, the thermal line- width of ^^C0(2-l) 
is 0. 15 km s"^ , which closely agree with the observed median ve- 
locity dispersion of this nearby cloud. Hence, Rinear resembles 
a cold and mainly starless cloud with an almost thermal velocity 
dispersion. 

3.2.3. The Coalsack cloud R2 comprising Tapia's Globule 2 

Figure [TOl presents the 1st and 2nd moment maps (intensity- 
weighted peak- velocities and velocity dispersion) for the south- 
weste r n region R2 , also k nown as "Tapia' s Globu le 2" (iTapiaL 
119731: iLada et al.L 120041: JRathbome et aD. l2009h . The peak- 
velocity distribution in this region does not exhibit as clear a 
velocity gradients as seen in Rlfar and Rinear (see previous sec- 
tions), but we find velocity gradients from the east and west of 
the main globule toward the extinction peak. East and west of 
the main peak, the 1st moment map reveals peak velocities — 
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5.8 km s"^ whereas in a relatively broad north-south filament 
extending almost across the central extinction peak, the peak 
velocities are found to be around -5.5 km s"^. Hence, in these 
^^C0(2-l) data we do not recover the east- west velocity gradi- 
ent previously discussed by lRathborne et al.l (,2009.) for the inner 
6' X 6' around the main extinction peak. 

The velocity dispersion of this region is also very narrow, 
spanning a range mainly below 1 km s"^ (Fig. [TO] bottom panel). 
In contrast to the Rl region where the average and mean values 
of the velocity dispersions deviated from each other, which indi- 
catis a broader velocity spread, for R2, the average and mean ve- 
locity dispersion are approximately the same with ~0.4kms"^. 
This can be interpreted as a sign of a relatively uniform velocity 
dispersion over the region. While the velocity dispersion at the 
edge of the map is around 0.2 km s"\ in the outskirts of the main 
extinction peak it rises to values on the order of ~0.4kms"^ 
and toward the extinction peak one finds values on the order of 
of 0.6 km s"^ or a bit higher. While even the latter values are 
still narrow compared to more active star- forming regions (e.g., 
[Goldsmith et al. 2008), they nevertheless are indicative of non- 
thermal gas motions in the direction of the highest extinction. 
For a more detailed discussion about the dynamical properties 
of the region, see section [42| 
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Fig. 10. Compilation of ^^C0(2-l) 1st and 2nd moment maps 
toward the southwestern region R2 in the top and bottom panel, 
respectively. The thick contours outline the extinction from 
iKainulainen et al.l (l2009h in 1 mag steps starting at an Ay of 
2 mags. Linear scale-bars are shown in the bottom-right corner 
of each panel. 



4. Discussion 

4.1. Fragmenting filaments 

The overall extent of the two clouds - "Nessie" discussed in 
I I I I 

Jackson et al. (2010) and Rlfar presented here - is similar with 
~80 pc and ~88 pc, respectively, and they are at similar distances 
(both at -3.1 kpc). However, there are also considerable diff'er- 
ences between the two clouds. First of all, the length-to-width ra- 
tio is diff'erent because while "Nessie" is extremely narrow with 
an approximate cylinder radius rcyi of only -0.5 pc, the structure 
within Rlfar is less well defined and hence it is also more dif- 
ficult to accurately determine a radius. However, based on the 
extinction maps, we roughly estimate an approximate average 
radius rcyi of ~ 2 - 2.5 pc for Rlfar, about four to five times 
larger than that of "Nessie", der ived from thei r MIPSGAL mid- 
infrared emission. Furthermore J Jackson et al.l (.2010.) report that 
star- formation activity is clearly identified toward many of their 
clumps, while our analysis reveals barely any star formation ac- 
tivity toward Rlfar. In addition to this, as discussed in section 
13.11 because the mass of Rlfar is distributed over such a large 
spatial extent, it cannot be considered as a "typical" very young 
high-mass star- forming region or IRDC. 

How does this structure compare to a gravitationally bound 
isothermal gaseous cylinder? Following Jackson e t al.l (12010.) a s 



well as the origi nal work by [C handrasekhar & Fermi[ ([19531) . 
[Nagasawa ('1987') and Inutsuka & Miyama (1992), we try to es- 
timate characteristic fragmentation scales for diff'erent condi- 
tions within self- gravitating gas cylinders. For an incompress- 
ible fluid, a characteristic length scale can be defined as /Ifrag ~ 
llrcyi, whereas for an infinite isothermal gas cylinder, the rela- 
tion is /Ifrag ~ 22// where H is the isothermal scale height given 
by // = Cs(47iGp)~^^^ with Cs the sound speed, G the gravita- 
tional constant and p the central gas mass density. 

Depending on the sound speed Cs we can calculate the 
characteristic scale height of such a filament. With a thermal 
sound speed of the gas at 15 K of Cs ~ 0.23 km s"^ and using 
^crit ~ 9x10^ cm"^ as a proxy for the density, the character- 
istic scale-height is // ^ 0.07 pc. However, inspecting Fig. [51 
we find that the region rather exhibits an average velocity dis- 
persion of ~1.26kms"\ which we use as a proxy for the full- 
width-half-maximum FWHM (see Sec. 13.2.11) . If one now re- 
places the sound speed Cs with the Gaussian velocity disper- 
sion cr = 1/2.35 ^ FWHM, we obtain a higher scale-height of 
H ^ 0.16 pc. This implies that the approximate cylinder radius 
Tcyi ~ 2.25 pc in any case exceeds the characteristic scale height 
by more than an order of magnitude. Because the characteristic 
length scale /Ifra g ^ H^cvi for an incompressible fluid requires 
Tcyi ^ H (^e.g. J Jackson et al .11201 Oh . the incompressible case can 
safely be excluded for Rlfar. 

Looking at the characteristic fragmentation scale for the in- 
finite isothermal gas cylinder /Ifrag ~ 22//, we derive /Ifrag ~ 
1.46pc or /Ifrag ~ 3.5 pc for the sound speed Cs or the Gaussian 
velocity dispersion cr, respectively. Because our data are clearly 
not described by a pure thermal line- width (Fig. [9]), the frag- 
mentation scale /Ifrag ~ 3.5 pc for a Gaussian velocity dispersion 
consistent with our data appears the most reasonable character- 
istic length scale we can derive in the framework of this picture 
of a gravitationally bound isothermal gaseous cylinder. 

Comparing this latter /Ifrag ~ 3.5 pc with our derived me- 
dian clump separation of ~4pc (Sec. 13.2. T1) , the data are broadly 
consistent with a compressible gravitationally bound isother- 
mal gaseous cylinder, in particular if one considers the observed 
spread in velocity dispersion values as well as an intrinsic spread 
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in the gas densities. However, they are inconsistent with an in- 
compressible fluid. 

Another way to characterize self-gravitating cylinders is in 
the framework of maximum critical ni a ss per u nit length M//max 
along the cylinder's axis (lOstrikej. 119641: iFiege & Pudritzl 
l200 Qa.b). Above the maximum critical linear mass density, 
the cylinders would be doomed to collapse. Following again 
iJackson et al.l (120101) and using the mean ^^C0(2-l) line- width 
of Av ^ 1.26 km s"\ for turbulently supported filaments we find 
a maximum critical M/ /max = 84(Av)^M0pc"^ = 133Mopc"^ 
(about a factor 10 higher than the maximum cr itical M/Zmax 7 
15Mopc"^ for undisturbed thermal filaments, Ostriker 1964; 
ilnutsuka & Miyama 1997 ). To compare that with our observa- 
tions, we can obtain a rough estimate of the observed Mjl by 
dividing the total mass estimated from the ^^CO data in Sec. 13. II 
(-2600 Mq) by the approximate length of the filament of ~73 pc 
in the ^^CO data. This results in an approximate observed Mjl 
of only ~36Mopc"^. Although this is larger than the maxi- 
mum critical linear mass density for thermal filaments (and also 
the measured Mjl ratio in, e.g., a thermal filament in Taurus, 
ISchmalzlet al. 2010), the observed mass-per-length scale is ap- 
proximately a factor 4 below the maximum critical M/Z^ax ratio 
for turbulently supported filaments based on the ^^C0(2-l) line- 
width. While a lower than critical M//max does not necessarily 
imply stability against collapse, it is at least consistent with the 
Coalsack at large not (or only barely) forming stars at the mo- 
ment. 

Comparing our results with those from the literature, we find 
that in the framework of isothermal gravitationally bound cylin- 
ders, we can distinguish for Rlfar between a compressible and 
an incompressible c onfiguration, wherea s this is not possible for 
the "Nessie" IRDC lJackson et al.1 (12010 ). Furthermore, the Mjl 
ratio for Rlfar with ~36Mopc"^ is considerably lower than that 
for ''N essie" or Orion (-l lOMppc"^ and ~385Mopc"\ respec- 
tively, Jackson et al. 2010). While one cannot directly compare 
the two latter values with the threshold we estimated above be- 
cause the observed line- width in these clouds are also larger than 
in the Coalsack, nevertheless, this comparison confirms the ex- 
pected evolution that starless clouds have low Mjl ratios that 
increase with time during the collapse and star-formation pro- 
cesses. However, not all low M// clouds have to evolve like that. 
They could also be transient structures that are dispersing again. 



4.2. Dynamical properties of R2/Tapla's Globule 2 



iRathborne et al.l 12009") recently suggested that the ring-like 
structure within this globule may have formed from the me rging 
of two sub-sonic flows. In contrast to lRathborne et al.l ("2009), we 
do not find a velocity gradient across the globule. This discrep- 
ancy likely arises from diff'erent obser vational biases, e . g., we 
observed the region in ^^CO, whereas IRathborne et al.l (l2009l) 
employed the rarer isotopologue C^^O. Furthermore, their C^^O 
map covered a much smaller area of the globule, and if one fo- 
cuses on the central few arcminutes of our data (Fig. \T0\ top 
panel), one could also get the impression of a velocity gradient. 
Nevertheless, the velocity dispersion increase toward the center 
is suggestive of dynamical evolution within the core. To further 
investigate the dynamics of the region, we produced a central- 
velocity-increment (CVI, e.g.. iHilv-Blant et al.l [2QQ8b map of 
the region, which is a diff'erent representation of the velocity 
changes over the observed region. Because the velocity structure 
is more prominent in east- west direction, the CVI map presented 
in Figure \TT\ studies the velocity changes in the R. A. direction 
with 10 pixel step sizes, corresponding to ~ 143'' or approxi- 



mately five resolution elements (section [2l). The observed peak- 
velocity changes are small toward the central peak of the extinc- 
tion map, however, we see four enhancements of the CVI-values 
toward the eastern and western edges in the north and south of 
the main extinction peak. This is consistent with the fact that we 
see a filamentary structure going from north to south through 
Globule 2 in the 1st moment map (Fig. [TOl) . Returning to the 
merging subsonic flow picture by Rathborn e et al. (2009), our 
data do not completely confirm this picture be cause we do not 
see a c lear east- west velocity gradient across as IRathborne et al.l 
(I2OO9I) do. However, the combined facts that we find clear CVI 
increases at the Globule's edges as well as an increase of veloc- 
ity dispersion toward the extinction peak are strongly suggestive 
of a dynamically active state of this dark globule within the large 
Coalsack complex. From these data we cannot establish the ex- 
act reason for the dynamical evolution, however, this Globule is 
an excellent candidate region for future active collapse and star- 
formation processes. 
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Fig. 11. ^^C0(2-l) central velocity increment (CVI) in R.A. to- 
ward the south western region R2. The th ick contours outline the 
extinction from lKainulainen et al.l (l2009l) in 1 mag steps starting 
at an Ay of 2 mags. A linear scale-bar is shown in the bottom- 
right corner. 



5. Conclusions 

Mapping approximately 1 square degree around the regions of 
the highest extinction in the almost entirely starless Coalsack 
dark cloud in ^^C0(2-l) with APEX revealed several interesting 
properties: 

- The Coalsack is not one coherent dark low-mass cloud at 
distances below 200 pc and typical velocities of ~ -4 to 
Okms"^ but we find a second velocity component at ~ 
-30kms"^ corresponding to a second cloud at a distance 
of ~ 3.1 kpc. This second more distant cloud dominates the 
dust-extinction signatures in the northeast of the Coalsack. 

- Mid-infrared MSX data reveal barely any star- formation ac- 
tivity toward any of the mapped clouds. 

- Although the total observed mass of the distant cloud with ~ 
2600 M© is far larger than those of the two mapped nearby 
clouds (on the order of several solar masses), we cannot con- 
sider it as a potential high-mass star- forming region, because 
its mass is distributed over a spatial extent in east- west direc- 
tion of ~73 pc. Hence, it is smoothly distributed over a large 
spatial area. 
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- The more distant cloud resembles in shape a large and Vassilev, V.,Meledin, D.,Lapkin, L, etal. 2008, A&A, 490, 1157 
twisted filament. Analyzing this structure within the frame- Williams, J. P., de Geus, E. J., & Blitz, L. 1994, ApJ, 428, 693 
work of an isothermal gravitationally bound cylinder, we 

find that the data are consistent with a compressible gaseous 
filament and inconsistent with an incompressible fluid. 
Furthermore, the characteristic mass per length M// is low 
with ~36Mopc"^. Hence, this structure is potentially stable 
against gravitational collapse. 

- The nearby low-mass clouds both exhibit narrow velocity 
distributions with median values between 0.2 and 0.4 km s"^ , 
also indicative of almost no star- formation activity. 

- Only Tapia's Globule 2 appears difl'erent. While the peak 
velocity dispersion barely exceeds 0.6kms"\ we neverthe- 
less see a line-width increase from the map edge toward 
the extinction peak. Furthermore, an analysis of the central- 
velocity-increments reveals significant velocity changes at 
the Globule's edges. Both features are suggestive of dynam- 
ical action in that region, e.g., maybe early infall activity. 

In summary, our data reveal diff'erent cloud components 
nearby as well as far away. While all observed clouds are mostly 
bare of any star- formation activity, we find filamentary structures 
consistent with compressible self-gravitating gaseous cylinders 
as well as one Globule in an apparent elevated dynamical state. 

Acknowledgements. This research made use of data products from the 
Midcourse Space Experiment. Processing of the data was funded by the Ballistic 
Missile Defense Organization with additional support from NASA Office of 
Space Science. This research has also made use of the NASA/IPAC Infrared 
Science Archive, which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National Aeronautics and Space 
Administration. 

References 

Beltran, M. T, Brand, J., Cesaroni, R., et al. 2006, A&A, 447, 221 

Brand, J. & Blitz, L. 1993, A&A, 275, 67 

Chandrasekhar, S. & Fermi, E. 1953, ApJ, 118, 116 

Churchwell, E., Babler, B. L., Meade, M. R., et al. 2009, PASP, 121, 213 

Corradi, W. J. B., Franco, G. A. R, & Knude, J. 1997, A&A, 326, 1215 

Fiege, J. D. & Pudritz, R. E. 2000a, MNRAS, 311, 85 

Fiege, J. D. & Pudritz, R. E. 2000b, MNRAS, 311, 105 

Fontani, F, Beltran, M. T, Brand, J., et al. 2005, A&A, 432, 921 

Franco, G. A. R 1995, A&AS, 114, 105 

Goldsmith, R R, Heyer, M., Narayanan, G., et al. 2008, ApJ, 680, 428 

Glisten, R., Nyman, L. A., Schilke, R, et al. 2006, A&A, 454, L13 

Hily-Blant, R, Falgarone, E., & Pety, J. 2008, A&A, 481, 367 

Inutsuka, S. & Miyama, S. M. 1992, ApJ, 388, 392 

Inutsuka, S. & Miyama, S. M. 1997, ApJ, 480, 681 

Jackson, J. M., Finn, S. C., Chambers, E. T, Rathborne, J. M., & Simon, R. 2010, 
ApJ,719, L185 

Kainulainen, J., Beuther, H., Banerjee, R., Federrath, C, & Henning, T 2011, 
A&A, 530, A64+ 

Kainulainen, J., Beuther, H., Henning, T, & Plume, R. 2009, A&A, 508, L35 

Kato, S., Mizuno, N., Asayama, S., et al. 1999, PASJ, 51, 883 

Klein, B., Philipp, S. D., Kramer, I., et al. 2006, A&A, 454, L29 

Knude, J. & Hog, E. 1998, A&A, 338, 897 

Lada, C. J., Huard, T L., Crews, L. J., & Alves, J. R 2004, ApJ, 610, 303 

Miettinen, O. & Harju, J. 2010, A&A, 520, A102+ 

Nagasawa, M. 1987, Progress of Theoretical Physics, 77, 635 

Nyman, L., Bronfman, L., & Thaddeus, R 1989, A&A, 216, 185 

Nyman, L. A. 2008, Handbook of Star Forming Regions, Volume 2 

Ostriker, J. 1964, ApJ, 140, 1529 

Price, S. D. 1995, Space Sci. Rev., 74, 81 

Rathborne, J. M., Lada, C. J., Walsh, W., Saul, M., & Butner, H. M. 2009, ApJ, 
690, 1659 

Reid, M. J., Menten, K. M., Zheng, X. W., et al. 2009, ApJ, 700, 137 

Rohlfs, K. & Wilson, T L. 2006, Tools of radio astronomy (Tools of radio as- 
tronomy, 4th rev. and enl. ed., by K. Rohlfs and TL. Wilson. Berlin: Springer, 
2006) 

Schmalzl, M., Kainulainen, J., Quanz, S. R, et al. 2010, ApJ, 725, 1327 

Tapia, S. 1973, in lAU Symposium, Vol. 52, Interstellar Dust and Related Topics, 
ed. J. M. Greenberg & H. C. van de Hulst, 43 



